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Introduction 

ON a swept wing, contamination along the leading 
edge, Tollmien-Schlichting waves, stationary or traveling 
crossflow vortices, and/or Ta.ylor-G6rt.ler vortices can 
cause the catastrophic breakdown of laminar to turbu- 
lent flow, which leads to increased skin-friction drag for 
the aircraft. The discussion in this Note will be limited 
to disturbances which evolve along the attachment line 
(leading edge of swept wing). If the Reynolds number of 
the attachment-line boundary layer is greater than some 
critical value, then the complete wing is inevitably en- 
gulfed in turbulent flow. Essentially, there are two crit- 
ical Reynolds number points that must be considered. 
The first, is for small-a.mplitude disturbances, and the 
second is for bypass transition. 

Summarized in Table 1, the experimental and the- 
oretical results agree for the critical Reynolds num- 
ber where small-a.mplitude disturbances become unsta- 
ble on the attachment, line. 1-4 Accounting for all lin- 
ear terms, and using an eigenvalue problem approach, 
Hall et a.l. 5 studied the linear stability of disturbances 
in the attachment-line boundary-layer flow called swept. 
Hiemenz flow, which is sketched in Fig. 1. By assuming 
instability modes which were periodic along the attach- 
ment. line, the calculations by Hall et. a.l. 5 agreed with the 
experiments and with the direct, numerical simulations 
of Spalart b , Theofilis', Jimenez et. a.l. 8 , and Joslin. 9,10 

For la.rge-a.mplit.ude disturbances, turbulence de- 
cays below some critical Reynolds number and transi- 
tion to turbulence will occur above this point.. At. this 
critical point., termed bypass Reynolds number, transi- 
tion bypasses the conventional linear instability break- 
down process. Summarized in Table 2, the experiments 
show that disturbances are damped for Rg < 100 and 
the flow becomes turbulent, for Rg > 100. (Refs. 1 and 
11-14). 

Hall and Malik 15 attempted to explain this discrep- 
ancy between linear theory and t.lnl turbulent, suppres- 
sion limits by studying the nonlinear disturbances using 
weakly nonlinear theory and temporal DNS. Subcriti- 
ca.l instability was observed in the computations; how- 
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ever, this sub critical growth did not. provide the connec- 
tion between linear instability and the contamination 
regions. 

Note, the wide gap between the linear critical 
Reynolds number of Rg ~ 245 and the turbulent, sup- 
pression critical Reynolds number of Rg ~ 100. Bridg- 
ing this gap is important, for wing design. The present, 
study will use direct, numerical simulations to validate a. 
linear 2D-eigenva.lue prediction method based on parab- 
olized stability equations by Lin and Malik 15 . This 
method is considered because it. suggests that a. num- 
ber of symmetric and asymmetric modes exist, and are 
stable or unstable on the attachment, line depending on 
the Reynolds number. If validated, the approach would 
predict, a. number of modes which are linearly damped in 
the Reynolds number regime 100 to 245; however, these 
modes may grow nonlinea.rly and provide an explanation 
to this region. 


Table 1. Critical Reynolds numbers 
for attachment -line instabilities. 


Experiment. 

Critical Rg 

Cumpst.y A Head 1 

245 

Pfenninger A Bacon 2 

240 

Poll 3 - 4 

230 

Calculations 

245 


Table 2. Experimental critical points for 
attachment-line turbulence suppression. 


Experiment. 

Bypass Rg 

Pfenninger 11 

100 

Gregory A Love 12 

95-98 

Caster 13 

88-104 

Cumpst.y A Head 1 

100 

Poll 14 

100 


u 

\ 



Fig. 1. Sketch of attachment -line region 
of swept Hiemenz flow. 


1 




AIAA J., vol. 34, no. 11, 1996, pp. 2432-2434 


Problem Formulation 

In general, the velocities u = (u,v,w) and the 
pressure p are solutions of the incompressible, unsteady 
Navier-Stokes equations. The instantaneous velocities 
u and the pressure p may be decomposed into base-flow 
velocities U_ = ( U , V, W) and pressure P and disturbance 
velocities u = ( u,v,w ) and pressure p. The base flow, 
originally described by Hall et al. 5 , is referred to as a 
swept Hiemenz flow. The disturbance is found by solv- 
ing the 3D incompressible Navier-Stokes equations (in 
disturbance form) 

c)u 1 

-= + (M-V)«T(iZ-V)MT(M-V)i7 = — Vp+ — V 2 m (1) 

dt R 

V • u = 0 

with boundary conditions 

u = 0 at y = 0 and u —>■ 0 as y —>■ oo (2) 

A Cartesian coordinate system x = ( x,y,z ) is used in 
which x is aligned with the attachment line, y is wall nor- 
mal, and z corresponds to the direction of flow acceler- 
ation away from the attachment line. A boundary-layer 
thickness is defined in the t/z-plane as 8 = \J vL/W 0 ; a 
Reynolds number, as R = IJ 0 b /v, U 0 ,V 0 , W 0 are velocity 
scales, and L is the length scale in the flow-acceleration 
direction z. If the attachment line is assumed to be in- 
finitely long, the velocities become functions of z and y 
only, and the similarity solution can be found. 

Numerical Method of Solution 

The 3D DNS code as described by Joslin 10 is used 
for the present validation study. High-order finite- and 
central-differences are used in the attachment-line (x) 
direction, and Chebyshev series are used in the wall- 
normal (y) and flow-acceleration (z) directions. For time 
marching, a time-splitting procedure was used with im- 
plicit Crank-Nicolson differencing for normal diffusion 
terms; an explicit three-stage Runge-Kutta method was 
used for the remaining terms. 

The disturbances are assumed to be from the dis- 
crete spectrum, which exponentially decay with distance 
from the wall. At the wall, in the far Held, and at 
the flow-acceleration boundaries, homogeneous Dirich- 
let conditions are imposed. The base flow is used for 
the inflow boundary condition. 

Symmetric and Asymmetric Instability Approach 

Recently, Lin and Malik 16,17 has shown with the- 
ory that both symmetric and asymmetric instabilities 
are present in incompressible and compressible swept 
Hiemenz flow. The solutions posed by Lin and Malik 16 
took the form 

{u,v,w}(x,y,z,t) = {u,v,w}(y,z)e l(ax -^ (3) 


Substituting this form into the Navier-Stokes equations 
leads to a system of partial differential equations in 
the flow-acceleration and wall-normal directions. The 
boundary conditions for the z boundaries took the form: 
Symmetric: 


du 

dz 


dv 

dz 


= w = 0 


at z = 0 


(4) 


{u,v}(y,z) = {u,v}(y,-z), w(y,z) 


at 


Z — Zmax 


w(y, -z) 


(5) 


Asymmetric: 


div 

u = v = — — = 0 
oz 


at z = 0 


( 6 ) 


{u, v}(y, z) = -{u, v}(y, -z), w(y, z) = w(y, -z) 

at z = z max (7) 

For the simulations, the entire attachment-line region is 
included within the computational domain, and there- 
fore, the boundary conditions at z = 0 are not needed. 

The theory suggests that the most unstable modes 
follow the sequence: symmetric (SI), asymmetric (Al), 
symmetric (S2), etc. where the growth rates of modes 
are 51 > Al > S2 > A2 > S3..., without exception. 
This theory and modal growth ordering were recently 
confirmed by Fedorov 18 using an asymptotic theory. Al- 
though according to the Fedorov analysis, the validation 
of a single mode implies the validation of all modes, here 
the first two dominant modes are simulated. 

From the results of the Lin-Malik technique, the 
wavenumber and growth rate for the first three modes 
at R = 700 (Rg ~ 282) and us = 0.1017 are shown 
in Table 3. The simulation of a pure mode will prove 
difficult because the discrimination of the wavenumbers 
would be difficult. The theoretical results suggest that 
the previous simulations of ’’discrete modes” are in fact 
spectrally rich. To use suction and blowing to generate 
the SI mode in the absence of the S2 mode is proba- 
bly not possible. However, a discriminating factor can 
be attributed to the phase relation between the sym- 
metric versus asymmetric modes across the attachment 
line and in the flow-acceleration direction. This differ- 
ence is obvious from the z = 0 boundary conditions (4) 
and (6). Hence, simulations could discriminate between 
symmetric and asymmetric modes. 
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Table 3. Lin and Malik lb eigenvalues 
for swept Hiemenz flow 
at R=700 and lo = 0.1017. 


Mode 

a r 

Ol i 

SI 

0.27481152 

-0.226959 x 10“ 2 

Al 

0.2751524.3 

-0.105988 x 10“ 2 

S2 

0.27548905 

±0.148157 x 10 -3 


Results 

The simulations are performed on a grid of 661 
points (~ 60 points per wavelength) along the attach- 
ment line, 81 points in the wall-normal direction, and 25 
points in the flow-acceleration direction. The far-held 
boundary is located at 40(5 from the wall, the computa- 
tional length along the attachment line is 216(5, and the 
flow-acceleration boundaries are located ±100(5. The to- 
tal Cray C-90 cost, for each simulation is 13 hrs for 8 pe- 
riods in time. Separately, the symmetric (SI) and asym- 
metric (Al) modes wert forced using suction and blow- 
ing. The phase of the Al mode in the flow-acceleration 
direction was determined using the Lin-Malik technique. 

In Figs. 2 and 3, the simulation results are com- 
pared with the wave growth rates described by the the- 
ory (listed in Table 3). The agreement is remarkably 
good considering the differences between the DNS and 
assumed solution form (4). For the theory, the Al 
mode has a constant wavenumber and growth rate in the 
flow-acceleration direction, while the simulations has & 
truly three-dimensional instability, and therefore, spec- 
tral differences in the r-direction are inevitable in this 
3D flow. To make this comparison, the results for the 
simulation are averaged over the flow-acceleration sta- 
tions: z = 0 and z = ±6. These stations were selected 
because:;,, as Joslin 10 shows, the streamlines very near 
the attachment-line are essentially aligned with the 2D 
attachment-line flow. The z = ±6 stations permit a 
cancelation of any opposing flow-acceleration effects. 

Concluding Remarks 

In this study, results are presented for the spatial di- 
rect numerical simulations (DNS) of 3D symmetric and 
asymmetric disturbances that propagate along the at- 
tachment. line of swept. Hiemenz flow. The comparison 
between the DNS results and the calculations of Lin and 
Malik 113 ’ 1 ' demonstrate that, both symmetric and asym- 
metric modes are present, in the attachment-line flow 
and that the theory adequately predicts these modes 
and the relative dominance of each mode. 

Although the connection between the linear in- 
stability Rg ~ 245 and the turbulent. Contamination 
Rg ~ 100 regions was not. definitively explained, it. is 
clear that. a. wealth of instabilities can be present, within 
this 3D flow field and that, some combination of these 


modes interacting in a. nonlinear manner will likely re- 
solve this regioaoof study. 
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Fig. 2. Symmetric disturbance (SI) growth 
compared with theory for three-dimensional 
attachment -line basic flow for R = 700 and u> = 
0.1017. 

10.0 


7.5 


vj 5.0 


2.5 
0.0 

0 50 100 150 200 250 

X 

Fig. 3. Asymmetric disturbance (Al) growth 
compared with theory for three-dimensional 
attachment -line basic flow for R = 700 and u> = 
0.1017. 
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